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Abstract 

Introduction. The article studies the problem of investigation of coupled nonstationary thermoelectroelastic fields in 


piezoceramic structures. The main approaches related to the construction of a general solution to the initial non-self- 
adjoint equations describing the process under consideration are briefly outlined. The work aims at constructing a new 
closed solution to the axisymmetric thermoelectroelasticity problem for a long piezoceramic cylinder. 

Materials and Methods. A \ong hollow cylinder whose electrodated surfaces were connected to a measuring device 
with large input resistance was considered. On the cylindrical surfaces of the plate, a time-varying temperature was 
given. The hyperbolic theory of Lord—Shulman thermo-electro-elasticity was used. The closed solution is constructed 
using a generalized method of finite integral transformations. 

Results. The developed calculation algorithm makes it possible to determine the stress—strain state of the cylinder, its 
temperature, and electric fields. In addition, it becomes possible to investigate the coupling of fields in a piezoceramic 
cylinder, as well as to analyze the effect of relaxation of the heat flow on the fields under consideration. 

Discussion and Conclusion. The use of assumptions about the equality of the components of the temperature stress 
tensor and the absence of temperature effect on the electric field allowed us to formulate a self-adjoint initial system of 


equations and construct a closed solution. 
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Introduction. Recently, different-purpose technical devices made of piezoceramic material have become 
widespread. Here, devices whose operation is based on the effect of the coupling of elastic, electric and temperature 
fields hold a specific place [1]. Various theories of thermoelectroelasticity were developed to describe their work taking 
into account the coupling of fields [2-4]. At the same time, for a better description and evaluation of non-stationary 
processes in structures, it was required to construct analytical solutions. However, the mathematical formulation of the 
problems under consideration included a system of non-self-adjoint partial differential equations, whose integration was 
difficult to treat mathematically. 

To solve this problem, as a rule, equations are investigated in an uncoupled form [5, 6], infinitely long bodies are 


analyzed [7—11], or thermoelectroelasticity problems are considered in a quasi-static formulation [12, 13]. 
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In this paper, we consider a coupled dynamic thermoelectroelasticity problem for an infinitely long hollow 
piezoceramic cylinder. As a result of the transformation of the initial calculated ratios, it was possible to form a self- 
adjoint system of equations, the integration of which was carried out by the method of incomplete separation of 
variables in the form of a generalized finite integral transformation [13]. 

Materials and Methods. Let a hollow, long, loose in the radial plane, piezoceramic cylinder occupy area © in the 


cylindrical coordinate system (r,. 0, z) : {a <r,<b,0<S0<21,-10<z< oo} . On the cylindrical surfaces, the 
temperature is given in the form of the following nonstationary functions (boundary conditions of the Ist kind) — 
@:(t.) (4, =a), 5(t.) A =); The internal electrodated surface is grounded, and the external one is connected to a 
measuring device with a large input resistance (electric idle mode. 

The mathematical formulation of the axisymmetric problem under consideration in dimensionless form includes 
differential equations of motion, electrostatics, and thermal balance based on the hyperbolic Lord—Shulman theory, as 


well as the boundary conditions [2, 7, 15]: 
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U'(r,,t.),0 (7,t.),O' (%,t.) — accordingly, the radial component of the displacement vector, the electric field 


potential, and the temperature increment in dimensional form; ( ©" (7,,t,)=T7(7,,t.)-T)(%), — current temperature 


as? P> mss €33— elastic moduli, density, piezo module, and 


t 


and temperature of the initial state of the body; C 
permittivity coefficient of electroelastic anisotropic material; (m,s = 3 V11>7%33— components of the temperature 


stress tensor (7,, =C,,@,,7%3; =C,,a@,); A,k,a,— coefficients of thermal conductivity, volumetric heat capacity, and 
linear thermal expansion of the material; g,— component of the pyroelectric coefficient tensor, £,,,— relaxation time; 


: ie 0 1 
©,— the rate of temperature change known at the initial moment; V = ap +-. 
r + 


In the case of grounding of the inner surface of a piezoceramic element, electrical voltage V(t.) is determined by 
the potential on its outer surface: 

V(t,)=¢(Lt.). (4) 

When constructing a general solution at the first stage of the study, the radial component of the electric field 


intensity vector is determined as a result of integrating the electrostatics equation: 
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U D, 
+d,—+a,0+—, (5) 
r r 


where D, — integration constant. 
Substituting (5) into (1)—(3) allows us to formulate a new problem with respect to functions U(r,t),@(r,t) .In 


this case, the condition of the absence of a radial component of the electric field induction vector on the outer 


cylindrical surface of the element (the last equality (2)) is fulfilled at D, =0; and the condition of grounding of 
the inner surface (Gr = 0) is satisfied as a result of integration (5). 


At the next stage of the solution, the inhomogeneous boundary conditions (2) are reduced to a form that allows 


further use of the procedure for incomplete separation of variables by the method of finite integral transformations. To 


do this, new functions u(r,t), N(r,t) related to U(r,t) ,O(r,t) are introduced: 
U(r,t)=A,(r,t)+u(r,t), O(7,t) =A, (r,t) +N(r,1), (6) 
where H,(r,t)=f,(r)A(LA+ 46 (r)A(R1)+ A (NO (+h (r)a, (0), 
H,(r.t)= f.(r)@(1)+ fo(")@, (1), 
f,(r)...fe(7) — twice differentiable function, A(r,t)_,, =(1+a,—a5—a,)U(r,1)/r- 


r=R,1 
Substituting (6) into the calculated ratios (1)-(3) with respect to functions U (r,t),O(7,1) when the following 
conditions are met: 


(1+a,)VH, +(a,—-1)H, = A(r,t) Hy,-2 =, Hy,. =, (7) 


r=R,1? 


allows us to get a new boundary value problem with respect to function u(r,t) ,N (r,t) : 
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It should be noted here that A(r,t) is a function of the displacements of the cylindrical surfaces of the cylinder. 


Initially, A(r,t) is equated to zero with its subsequent determination and refinement H,, F,, F,. 


Further transformations of the calculated ratios (8)—(10) are associated with the use of the following assumptions: 


b, =1,b, =0, and the introduction of a thermoelastic potential 


N=VB. (11) 
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Condition b,=1 can be accepted without a large error, since for piezoceramic materials b, =0.94+0.98, and 
dependence b, =0 is fulfilled in the case of equality of the components of the temperature stress tensor (7,, = 73, ) and 
the absence of temperature influence on the electric field ( g, =0). 


As a result, the following task is formed regarding u(r,t),B(r,t) ; 


1 ou 
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where F, = ae {4 pF Nan +a,H,),VH, = f;(r)@,(t)+ f(r) @ (t). 


The initial boundary value problem (12)—(14) is solved using the structural algorithm of the generalized finite 


integral transformation (FIT) [14]. At the same time, it is possible to use single-component unknown transformation 


kernel K (Ay r) for this task: 


{G,(/,,1),G,(A,,)} = [{u.), BU. )}KG, rar, (15) 


{u(r,t),B(r,t)} = iG (A,.t),G)(A,,.01K (Ar) Kl. (16) 


=1 


1 
|K, |’ = [KG,, r)’rdr ; 
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where A, — eigenvalues forming a countable set. 


As a result of using the FIT algorithm [14], we obtain problems with respect to the transformation kernel K (A, r) ; 


°K(A, IK (A, 
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The general solution to problem (17), (18) has the form: 
K(A 1) =% (A) A (Ar) —Jo (4 )¥ (Ar). (21) 
Here, eigenvalues A, are determined using the following transcendental equation: 
Y (A) Jo (A,R)-J, (4,)% (AR) =0. 
The system of differential equations (19) is reduced to the following resolving equation of the 4th order with respect 
to G(4,,0): 


d* d° d* 
+b 
= “dt dt? 


+s; + Dg, “ +b, |G, =F, (22) 
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Since the characteristic equation corresponding to (22), 
k* +b, +b k? +bk+b,, =0, 


is valid, then it, from the condition of the oscillating solution for G,(/,,f), has two real roots (k,,,k,,;) and two 


li? 
complex-conjugate roots 
(k,,=o0+i0, k,, =o-i@). 
In this case, the general solution to equation (22) has the form: 


G,(/,,t) = D,, exp(k,t) + D,, exp(k,,t) + D,, exp(k,t) + Ds, exp(k,,t)+ (23) 


wl Fa (r){exp[k, (t- t)| —exp| k,, (t —r) |} dt by, (z)exp| o(t 7) {Pro sin(wt—ar)—b,,,cos(at+ar)}dr, 
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Function G,(4,,f) is determined from the first equation of system (19). Substitution of the obtained 
expressions for the transforms under the boundary conditions (20) makes it possible to determine the integration 
constants D,,...D,;... (19). 


Substitution of G,(A,,t),G,(/,,f) into (16), (11), (6) allows us to get the final expressions for functions U (r,t) ‘ 


O(r,t): 


—2 


U(rnt)=H,(r.t)t 0G AK (ANIKI, (24) 


i=l 


O(r,t)= + La, VK (A,r) ||K, |. 


At the final stage of the study, functions H, (752), H, (7,1) are determined through solving the following differential 


equations: 


Fg epg SG, (25) 
or r or 


which makes it possible to significantly simplify the right parts (F, F,) of the calculated ratios (8). 


Mechanics 


85 


http://vestnik-donstu.ru 


86 


Advanced Engineering Research 2022. V. 22, no. 2. P. 81-90. ISSN 2687-1653 


Substitution of expressions for H,,H,in (25) enables to form systems of equations with respect to functions 


J, (r)...£,(r) , that are determined when the conditions are satisfied (7). 
The potential of the electric field of a piezoceramic cylinder is determined from integrating equality (5) and 


satisfying the next-to-last boundary condition (2): 


male OH, (r,t) +a, H, (r,t) +a ( rt) fr (26) 
or r 


* + De (t) 


K; 
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The obtained calculated relations (24), (26) satisfy differential equations (1) and boundary conditions (2), (3), i.e., 
they are a closed solution to the problem under consideration. 


Research Results. As an example, we considered a radially polarized piezoceramic cylinder (b =0.02m, R=0.8 ) 


of PZT—4 composition, having the following physical characteristics [10]: = 7500 kg/m* , 


{Ci1sCyssC,s} = (13.9; 11.5; 7.43} x10" Nim?, {e,,,€,}={-5.2; 15.1} Chm’, 
{Hiss Yeat = (4.6; 3.9} 10° H/(m* °C), €,, =5.62x10°’ ®/m, g, =2x10*Ka/(m’ °C), 


k =3x10° J m?* °C), A=1.6 W/(m °C), B,,, =5x10°c. 
A temperature load acts on the inner surface (7, =a ) of the piezoceramic cylinder: 


a, (Aly ae in( Et | -1,)+H(t, — tr.) Oo; (t,) =0, 


where H(?) —Heaviside step function (H(7)=1 at 720, H(f)=0 at 7 <0), 7, =Te —Tys Trhagsl nox — 
maximum value of the external temperature effect and the corresponding time in dimensional form 
Ge, 3373 K (100°C), 7, = 203-K (20°C), F.=1 ©), 

Figure 1 shows graphs of changes in functions ©* (a), U (24), (r,t) along radial coordinate r at various points 


in time t. The numbers 1—3 respectively indicate the results obtained at the following time values: t =1,,,.54t ax >! Snax 


Analysis of the calculation results allows us to draw the following conclusions: 

— sufficiently large value of the coefficient of linear thermal expansion a@, of the piezoceramic material causes rapid 
heating of the cylinder; 

— radial displacements on the inner cylindrical surface (r=R) at the first stage of the study (t=1,,,,) take the 
greatest values, followed by a decrease over time. The reverse pattern is observed with respect to the displacements at 


r=1; 


ry 
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@' (7,1), °C 


c) 
Fig. 1. Graphs of ©* (7,1), U(r,t),A(7,t) changes along radial coordinate r 


at various points in time t (1—t,4.52—Afinax >3—15tingx ! d— O" (r,t)+r 3 b— U(r,t)=r;c— O(r,t)+r 


The degree of coupling of thermoelectroelastic fields is most conveniently analyzed using the coefficient 


2 
by, = z( a,t+ 2 from equality (21). Here, a, determines the coupling of electroelastic fields, and “ _ the 
a, a, 
effect of the rate of change in the volume of the body on its temperature field. 
Figure 2 shows displacement graph U (L t) in time f taking into account (solid line) and without taking into 
account (dotted line) the induced electric field. 
It should be noted that the preliminary polarization of piezoceramics causes the formation of a more “rigid” material 


(a, = 0.353) and, accordingly, a decrease in displacements during deformation of the cylinder. 
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The coupling of temperature and electroelastic fields in a piezoceramic cylinder can be neglected due to the small 


value 


%s -1.8x107 «1. 
a, 


U(r,t) «10° 


0 1x10° 2x10° 3x10° 


Fig. 2. Graph of change U (1, t) in time ¢ 


(solid line — including the polarization; dotted line — without including the electric field) 


Figure 3 shows graphs of changes in electrical voltage V(t) over time, taking into account (solid line) and without 


taking into account (dotted line, £,, = 0) the heat flow relaxation. 


V(t) x 10° 


Fig. 3. Graph of change V(t) over time t(t,,, =0.001 s,t,,. =196) 


The calculation results show that for the problem under consideration, the refined hyperbolic Lord-Shulman theory 


a) >5.56x10° K/’and at 
dt. s 


|*=0 


should be used at a high rate of change in the temperature load ( t.,, > 0.001 s, 


lower speeds — the classical theory of thermoelectroelasticity ( £,,, =0). 

Discussion and Conclusion. The constructed new closed solution to the coupled dynamic problem with satisfaction 
of the boundary conditions of thermal conductivity of the Ist kind made it possible to determine all the components of 
thermoelectroelastic fields in a long piezoceramic cylinder. The advantage of the presented calculation algorithm is that 
there is no need to approximate the temperature function when studying the equation of motion, in contrast to the 
uncoupled formulation of the problem. At the same time, the effect of the rate of change in the volume of a 


piezoceramic body on its temperature field can be neglected. 
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